In-plane magnetic field anisotropy of the FFLO state in layered superconductors 
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There are strong experimental evidences of the Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) state 
formation in layered organic superconductors in parallel magnetic field. We study theoretically the 
interplay between the orbital effect and the FFLO modulation in this case and demonstrate that 
the in-plane critical field anisotropy drastically changes at the transition to the FFLO state. The 
very peculiar angular dependence of the superconducting onset temperature which is predicted may 
serve for unambiguous identification of the FFLO modulation. The obtained results permit us to 
suggest the modulated phase stabilization as the origin of the magnetic-field angle dependence of 
the onset of superconductivity experimentally observed in (TMTSF)2C104 organic conductors. 
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Layered superconductors exposed to an external mag- 
netic field aligned parallel to their conducting layers have 
been in the focus of theoretical and experimental investi- 
gations due to their remarkable anisotropic propertiesi*^ 
favorable to the formation of the spatially modulated 
phase, known as the Fulde-Ferrell-Larkin-Ovchinnikov 
(FFLO) state.— »^ In particular, in the family of organic 
layered superconductors (TMTSF)2X, where anion X is 
PFg, CIO4, etc., very large upper critical fields, which 
exceed the Pauli paramagnetic limit, were reported.—"— 
In layered conductors the orbital motion of electrons 
is mostly restricted to the conducting crystal planes 
when hopping between adjacent layers is small. Thus 
the magnetic field applied parallel to the conducting 
planes causes only small diamagnetic currents and the 
orbital depairing is strongly weakened. Therefore, spin- 
singlet superconductivity is mainly limited by the Zee- 
man energy (Pauli spin polarization) of the quasipar- 
ticlcs. In contrast, the Pauli effect is negligible for a 
spin-triplet pairing because in this case Cooper pairs gain 
Zeeman energy without loosing the condensation energy. 
The question concerning the singlet or triplet symme- 
try of the superconducting order parameter in layered 
organic conductors is a current topic of debate. In- 
deed, the Nuclear Magnetic Relaxation (NMR) exper- 
iments with (TMTSF)2PF6 salts below Tc and under 
pressure showed the absence of the Knight shift, thus 
supporting the triplet scenario of pairing^ii while the 
^''Se NMR Knight shift in a recent experiment with 
(TMTSF)2C104 revealed a decrease in spin susceptibil- 
ity consistent with singlet pairingJ^ ^■^C NMR measure- 
ments with K-(BEDT-TTF)2Cu(NCS)4 also evidenced 
for a Zeeman-driven transition within the superconduct- 
ing state and stabilization of FFLO phase For the sin- 
glet superconductivity the FFLO phase can be a candi- 
date for the enhancement of the upper critical fieldi^i^ 
Note that in the compound (TMTSF)2C104 the sub- 
stantial anisotropy within the conducting a — b* plane 
is present. When magnetic field is aligned along the 
high conductivity a-axis the orbital currents are strongly 
quenched, which favors the FFLO phase appearanceJ^ii^ 



Interestingly, for a magnetic field applied along the b*- 
axis, the 3D — > 2D dimensional crossover occurs in the 
high field regime and the coexistence of the hidden reen- 
trant and FFLO phases can emergeJ^ 

Recently the in-plane angular dependence of up- 
per critical field, Hc2, of the organic superconductor 
(TMTSF)2C104 has been measured for wide tempera- 
ture intervals The observed upturn of the Hc2 curve at 
low temperatures has often been discussed in connection 
with the possibility of the FFLO state formation] ^^1^^ In 
addition as shown in Ref. [l^ the superconducting phase 
in high magnetic field is more strongly suppressed by im- 
purities than that in low field, as expected in the FFLO 
scenario.— Furthermore, an unusual in-plane anisotropy 
of Hc2 in the high-field regime was observed, which was 
again interpreted as an evidence of the FFLO state sta- 
bilization. This argument is based on the prediction of a 
very peculiar in-plane angular dependence of the FFLO 
critical field due to the orbital effects in thin supercon- 
ducting filmsi^ Motivated by these experimental find- 
ings we investigate in this work the influence of the spa- 
tially modulated superconducting phase on the in-plane 
anisotropy of the upper critical field in layered supercon- 
ductors with s-wave pairing. 

To describe the layered superconductors we consider 
a system of layers in xy-plane, stacked along the z-axis. 
The single-electron spectrum is approximated by 
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where p = {px,Py,Pz) is the electron momentum. The 
in-plane motion is described within the effective mass 
approximation while the tight-binding approximation is 
used to describe the motion along the z-direction. The 
corrugation of the Fermi surface due to the coupling be- 
tween adjacent layers (interlayer distance d) is assumed 
to be small, i.e. t ^ Tco, but sufficiently large to make 
the mean-field treatment valid, \\n{Tco/t)\Tco/ Ep <C 
1.— Here TcO is the critical temperature of the sys- 
tem at zero magnetic field and Ep is the Fermi en- 
ergy. We choose a gauge for which the vector potential 
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A = H X r [r = (a;, y, 0) is a coordinate in xy-p\a.ne], i.e. 
Az = —xH sin a + yH cos a, where a is the angle between 
the applied field, with amplitude _ff , and x-axis. As was 
demonstrated in Ref. [13 the anisotropic model with effec- 
tive masses can be reduced to the isotropic one by a scal- 
ing transformation and corresponding renormalization of 
the magnetic field. Therefore, in the pure Pauli regime, 
the orientation of the FFLO modulation vector, q, is ar- 
bitrary in the case of an elliptical Fermi surface. However 
any deviation of the Fermi surface from the ellipticity 
fixes the direction of the modulationi^i^i^i Hereinafter 
we assume for the sake of simplicity that these devia- 
tions from ellipticity are small and their role is just to 
pin the direction of the vector q, which is supposed to 
be along the i-axis. Performing this scaling transforma- 
tion, we will thus consider from now on an isotropic in- 
plane spectrum, with mass m = nix, and a magnetic field 
H = iJ[(ma;/mjy)^/^ cosa, sina, 0]. Taking into account 
that the system is near the second-order phase transi- 
tion, the linearized Eilenberger equation on the anoma- 
lous Green function /a;(n, r^p^) describing layered super- 
conducting systems acquires the form (for positive Mat- 
subara frequency uj at temperature T)^ 

[w+i/i+ivF.V+2it sin(p^d) sin(Q.r)]/^(n, r,^^) = A(r). 

(2) 

Here h = hbH is the Zeeman energy, xp = 
vp'^ is the in-plane Fermi velocity, and Q = 
(7r(i7?/(?!)o)[~ sina, (maj/my)-'-/^ cosa, 0] with 0o ~ 7rc/e. 
The order parameter is defined self-consistently as 
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where A is the BCS pairing constant and the brackets 
denote averaging over pz and n. Here we considered a 
layered superconductor in the clean limit, meaning that 
the in-plane mean free path is much larger than the corre- 
sponding coherence length, = ^^F/(27rTco). The upper 
critical field corresponds to the values of H for which the 
system of Eqs. ([2]) and ([3]) can be solved. 

The solution of the Eilenberger equation © can be 
chosen without loss of generality as a Bloch function 
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Equation (U) takes into account the possibility for the 
formation of the pairing state (k-|- ^,t;— k-|- with 
finite center-of-mass momentum. At the same time, the 
order parameter can be expanded as 
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It is knowr^ that in the absence of orbital effect, the 
FFLO state only appears at T < T* ~ 0.56Tco and H > 
H* ~ 1.06Tco///B, where (T*, H*) is the tricritical point. 
Therefore, the order of the magnitude of the magnetic 



field required to observe the FFLO state can be found 
from the relation fisH ^ T^o. Taking this into account 
one obtains vpQ ^ vpedTca/ jjlbC ~ {d/a)Tco, where a is 
the unit cell in the xy -plane. Therefore vpQ > Tcq. Due 
to the assumption t <C Tco ^ vpQ one has \/tTco ^ vpQ. 
This condition allows us to retain only the terms up to the 
first harmonics in Eqs. (U) and ([S]) [since we will retain 
only the terms up to (t/Tco) in the final expressions]. 
Substituting Eqs. (g]) and ([5]) into Eq. ^ one gets 



i(q)/o + <7-i = Ao, 

i(q±Q)/±i±i7o = 0, 
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where L{q) = u + ih + ivp.q/2 and t = tsin{pzd). If 
one neglects the Zeeman term these equations readily 
describe the reentrant phase predicted by Lebed;^ with 
critical temperature Tco at fields H ^ (l>o/{d£,o)- While 
keeping the terms up to the second harmonics within 
the same procedure would yield the Lawrence-Doniach 
equationi^ Inserting the solution of Eqs. ([H]) into the self- 
consistency equation ([3]), keeping only terms up to the 
second order in t/Tco, and subtracting it with a similar 
equation relating A with Tco, we obtain 
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L2(q)Vi(q + Q) i(q-Q) 

This equation defines the temperature dependence of the 
upper critical magnetic field H^2 i^i layered superconduc- 
tors, when both the paramagnetic and orbital effects are 
accounted for. 

In the limit i <C Tc the magnitude of the FFLO modu- 
lation vector can be calculated by neglecting the orbital 
part in Eq. ([7]). When making average over the Fermi 
surface one gets the equation 
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with reduced variables h = h/inT and q = q/2'nT, which 
gives rise to a FFLO vector q with magnitude that max- 
imizes the upper critical field, thus defining Tcp{H) and 
qp{H), in the pure Pauli limiti^ Finally averaging Eq. 
([7|) over the Fermi surface one obtains the equation for 
the onset of superconductivity in layered superconduc- 
tors, Tc{H), in the presence of both Zeeman and orbital 
effects: 
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The summation over the Matsubara frequencies is per- 
formed numerically. We used A'^ = 10^ terms in the 
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FIG. 1; Reduced critical temperature (with respect to the critical temperature in the pure paramagnetic limit) as a function 
of in-plane magnetic field in a layered superconductor, for several angles a between the field and i-axis (equivalently FFLO 
modulation vector at _ff > H*) Left panel: isotropic regime with rrix — rUy. [Inset: (_ff, r)-phase diagram in the pure 
paramagnetic limit.] Right panel: anisotropic regime. Solid lines are for nix = Wrriy; dashed lines are for tUx = O.lmy. The 
calculations are performed for Fermi velocity vp ~ 2.10^ cm.s~^,— and interlayer distance d = 1.3 nm.— 



summation and this number suffices for convergency at 
T/Tco > 10~^. Fig. [T] sfiows tlie variation of the nor- 
malized correction of the transition temperature, AT^ = 
Tc — Tcp, as a function of reduced strength of the mag- 
netic field, H/HpQ and angle a (between H and x-axis). 
Here Hpo = Ao//i_B is the critical magnetic field at T = 
in Pauli limited two-dimensional superconductors.— [The 
(H, r)-phase diagram in this regime is given in the inset 
of left panel.] The left panel describes the isotropic sit- 
uation, typical for layered quasi-2D compounds,— while 
the right panel exhibits results obtained for the highly 
anisotropic in-plane Fermi surface of layered conductors, 
exhibiting quasi-lD character.— We consider two oppo- 
site mass anisotropics. When rrix = lOmj, q is along the 
heavy mass direction, while in the case of — O.lmy 
it is along the light mass direction. As it was intu- 
itively expected the orbital effects reduce the supercon- 
ducting onset temperature, ATc < 0. While increasing 
the applied magnetic field, ATc first decreases in most 
cases until the tricritical point, H*, is reached and the 
curve of ATc exhibits a kink. At H > H* the func- 
tion ATc{H) strongly depends on the in-plane effective 
mass anisotropy and angle a. For a close to 90°, ATc 
exhibits an upturn and Tc approaches the paramagnetic 
limit, Tcp, when H increases. In contrast, for small a an 
increase of the magnetic field leads to a decrease of ATc. 
For intermediate angles, ATc can be a non-monotonic 
function of the field strength. In the isotropic case and 
for H/Hpo ^0.75 the largest correction to the onset tem- 
perature I ATc (o^)! occurs at a sa 20°. For mx/my = 10 
and H/Hpa > 0.8, |ATc'"'^''| is at angles close to a « 45°, 
while for nix /my = 0.1 and H/Hpo > H* , |ATc"'^^| is 
at angles close to a « 0°. One can infer that the strong 
field-direction dependence of the superconducting onset 



temperature Tc (a), appears at high magnetic fields when 
the FFLO state develops, while it is absent at low fields. 

The change in the anisotropy of the superconducting 
onset temperature that is induced by the FFLO phase is 
particularly visible in Figs. [2]and[31 where the magnetic 
field angular dependence of the normalized superconduct- 
ing transition temperature, Tc (a) /T^p, at constant mod- 
ulus of the in-plane magnetic field and t/Tco = 0.2, is 
plotted. In the polar plot the direction of each point 
seen from the origin corresponds to the magnetic-field 
direction and the distance from the origin corresponds to 
the normalized critical temperature, when the orbital de- 
structive effect is taken into account. We show here such 
dependence because this type of representation is essen- 
tially informative and was realized in the experiment 
For magnetic fields below H* and rrix = my one can 
see an expected isotropic behavior of the upper critical 
field. When increasing H above H* , a strong in-plane 
anisotropy of Hc2 develops, which remains and becomes 
essentially pronounced at high fields. In particular, rel- 
atively strong dips at a = ±18° and a = ±162° with 
small peaks at 0° and 180° develop with external mag- 
netic field for the case of the isotropic in-plane Fermi 
surface. The maximum transition temperature is for the 
magnetic field orientation perpendicular to the direction 
of the FFLO modulation vector along 2;-axis. 

The peculiar features of the anisotropy of the super- 
conducting transition onset with the orientation of the 
in-plane magnetic field strongly depend on the anisotropy 
of the Fermi surface. For example, comparing the curves 
in Fig. [21 one can notice the totally different in-plane 
field-direction dependence of Hc2 above the tricritical 
point when the largest electron mass is along the x- (solid 
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FIG. 2: Normalized superconducting transition temperature, 
Tc (a) /Tcp as a function of the angle between the directions of 
the applied magnetic field and the vector q for several values 
of H/Hpo and for nix ~ rUy. Here t/Tco — 0.2. 




FIG. 3: The same as in Fig. [2]but for mx = lOrriy (solid line) 
rrix — O.lmy (dashed line). 

line) or ^/-direction (dashed line). For an in-plane mass 
anisotropy nix = lOnij,, the positions of the dips are at 
angles a = ±50° and ±130°, respectively. The positions 
of the small peaks are the same as in the isotropic case, 
however they are strongly broadened. For opposite mass 
anisotropy mx/my = 0.1, the dips are strengthened and 



are found now at a = ±5° and ±175°. The fine structure 
at a = 0°,180° is relatively sharpened. Therefore the 
shape of the observed Tc(a) dependence could permit to 
determine the orientation of the FFLO modulation vec- 
tor. However the common feature of the field evolution 
of the in-plane upper critical field anisotropy observed 
in all considered cases is that the anisotropy becomes 
more pronounced with the field strength. Furthermore, 
the maximum critical temperature in the FFLO regime 
always corresponds to a = ±90°, that is H perpendicu- 
lar to q, irrespective to the effective mass ratio. This is 
in contrast with the usual behavior near Tco, where the 
critical temperature is maximal for the magnetic field 
oriented along the lightest mass axis. 

In conclusion, our results show that the FFLO mod- 
ulation strongly interferes with the orbital effect and 
provides the main source of the in-plane critical field 
anisotropy. The superconducting onset temperature is 
maximal for the field oriented perpendicular to the FFLO 
modulation vector. The change of the anisotropy of the 
critical field as well as of its fine structure may give 
important information about the FFLO state and un- 
ambiguously prove its existence. Our calculations sup- 
port the interpretation of the experimentally observed 
in-plane anisotropy of the onset of superconductivity in 
(TMTSF)2C104 samples as a realization of the FFLO 
state with the modulation vector close to b*-axis»^ How- 
ever, the compound (TMTSF)2C104 is in fact in the 
regime t > T^o [t - 2-7 K and T^q = 1.45 In 
this case the FFLO vector can be changed by the or- 
bital effect. Its orientation will result from the inter- 
play of the Fermi surface non-ellipticity, which favors 
pinning of q in certain direction, and the orbital effect, 
which prefers to orient q perpendicular to H. Neverthe- 
less we expect that the obtained results will be qualita- 
tively applicable in this case as well. We suggest that 
the predicted in-plane anisotropy of Hc2 can be observed 
in experiments with k-(BEDT-TTF)2Cu(NCS)2 organic 
superconductorsj^i*^ For this salt the angle-dependent 
magnetoresistance measurements^ provide the estimate 
of the interlayer transfer integral t k, 1 — 2 K, which 
is much smaller than Tco = 9.1 K, and the orbital ef- 
fect should only slightly change the modulation vector;^ 
In this work we have assumed s-wave superconductivity, 
however it is not an important ingredient in the present 
theory and we expect similar results in the case of li-wave 
pairing, which provides an additional source of pinning 
for the modulation vector.— 
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